Spectral properties and pseudo-gaps in a model with d-wave pairing symmetry 



Bumsoo Kyung 

Departement de physique and Centre de recherche en physique du solide. 
Universite de Sherbrooke, Sherbrooke, Quebec, Canada J1K 2R1 
(March 30, 2000) 

A model with d-wave pairing symmetry is studied by employing a non-perturbative sum rule 
approach. At low temperature the magnitude of a normal state pseudogap shows strong k or angle 
dependence well fitted by cos 2(j> form. With increasing temperature, the pseudogap closes at some 
critical angle 4> c and beyond this angle a single quasiparticle-like peak appears. The resulting Fermi 
surface is strongly temperature dependent. Both in the spectral function and the density of states, 
the pseudogap disappears in a manner that the spectral weight fills in the pseudogap instead of 
closing it with increasing temperature. All these features are qualitatively consistent with ARPES 
for underdoped cuprates. 

PACS numbers: 71.10.Fd, 71.27.+a 



The nature of an excitation gap in high-temperature 
superconductors has been one of the puzzling issues in the 
community of condensed matter physicists. Through sev- 
eral years of extensive experimental work, general consen- 
sus regarding the superconducting gap symmetry seems 
to be reached that the superconducting gap has mainly 
d-wave character with possibility of a small mixture of 
other angular momentum states, in contrast to con- 
ventional BCS superconductors with an isotropic s-wave 
gap. Recent discovery of a normal state pseudogap in un- 
derdoped cuprates has shed another side of the anoma- 
lous behaviors in the copper oxide superconductors. For 
these materials the low frequency spectral weight begins 
to be strongly suppressed below some characteristic tem- 
perature T* higher than T c . This behavior has been ob- 
served through various experimental probes such as pho- 
toemission, ^-pj specific heat, |J7[ tunneling, || NMR, 



M and optical conductivity. |10| In particular recent an- 
gle resolved photoemission spectroscopy (ARPES) 
and tunneling experiments Q indicate that the pseudo- 
gap phenomenon is closely related to pairing fluctuations. 
These measurements clearly exhibit that the normal state 
pseudogap has the same angular dependence and magni- 
tude as the superconducting gap and that often the only 
difference between the spectra in the pseudogap state and 
the superconducting state is in their linewidths. Typ- 
ically T* is much higher than T c and their doping de- 
pendence is qualitatively different. While T c decreases 
with underdoping, T* increases in contrast. This fea- 
ture suggests that T* does not follow T c characterized by 
long-range phase coherence, but instead some kind of a 
mean- field critical temperature Tmf- In spite of several 
possible scenarios such as the spinon pair formation with- 
out the Bose-Einstein condensation of holons, 1 1 - 13|] 
strong superconducting phase fluctuations [p^-^6[ and a 
magnetic scenario near the antiferromagnetic instability 
Jl7| , [i"8| and so on, at present there is no consensus in the 
origin of the pseudogap. 



For the past several years extensive theoretical ef- 
fort has been also made by several groups to under- 
stand this anomalous pseudogap behavior in the context 
of short range effective (attractive) interaction between 
electrons. This may be divided into two different classes 
of approach. In the first class, quantum Monte Carlo 
(QMC) simulations were made for the attractive Hub- 
bard model. |19-2(| Although this Hamiltonian is not 
a realistic model for understanding the complex physi- 
cal behaviors in the underdoped cuprates, it is believed 
to capture the important ingredient of the paring fluc- 
tuations in those materials. In the absence of a small 
parameter this numerical method has played an impor- 
tant role in understanding of the model, in spite of some 
uncertainties due to finite size effect and numerical an- 
alytical continuation. In this approach the one-particle 
spectral function as well as the density of states clearly 
show the precursor of the superconducting gap in the nor- 
mal state. In the second class, various assumptions and 
approximations are used. This class includes the effect of 
vortex phase fluctuations on the single-particle properties 
|l6f and the 'paring approximation' theory, J^J and T- 
matrix and self-consistent T-matrix approximations for 
the two-dimensional attractive Hubbard model. |2^-|3^j 
More recently the self-consistent T-matrix approach was 
also applied to a model with d x 2_ y i pairing. |34|,|35| The 
present approach is a variant of the T-matrix approxima- 
tion applied to a model with d-wave pairing symmetry. 

In a previous study, |36| a non-perturbative sum 
rule approach was developed for the attractive Hub- 
bard model by extending previous work on the repul- 
sive Hubbard model. fl37| It is found that in two dimen- 
sions, the mean-field transition temperature is replaced 
by a crossover temperature where the correlation length 
starts to grow exponentially. At sufficiently low temper- 
ature, a Kosterlitz-Thouless 0(2) transition should oc- 
cur, but it is not reproduced by that approach since it 
is in the 0(n = oo) universality class pj|. Nevertheless, 
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the agreement with Monte Carlo calculations is quanti- 
tative for both one and two-particle correlation functions 
over the whole range of parameters accessible by Monte 
Carlo calculations where it is found that in two dimen- 
sions, a pseudogap appears in the single-particle spectral 
weight p^j as well as in the density of states. ^(J Re- 
cent dimensional crossover study by Preosti et al. ]39| ] 
shows that the pseudogap effect is basically absent in 
three dimensions. In weak to intermediate coupling, the 
appearance of a pseudogap is traced back to the grow- 
ing critical pairing fluctuations in the low-temperature 
renormalized classical regime of the low-dimensional sys- 
tem. With increasing temperature, the spectral weight 
fills in the pseudogap instead of closing it. Furthermore, 
the pseudogap appears earlier in the density of states 
than in the spectral function. It was also noted that the 
qualitative features found in this study should apply to 
the d-wave case. In this paper we study in detail spectral 
properties and pseudogaps in a model Hamiltonian with 
d-wave pairing symmetry which is more appropriate for 
high-temperature superconductors. 

We consider on a two dimensional square lattice a 
simple model Hamiltonian which has a superconducting 
ground state with d-wave symmetry 



H = > etct cr 

fc,<7 
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where eg = —2t(cosk x + cosk y ) and N is the num- 
ber of lattice sites. We choose p as a d-wave sep- 
arable potential given as VH(k)H(k') where II(fc) = 

k 2 -k 2 

cosk x — cosky, or as 1.5264 cos20 = 1.5264 fc 2 +fc 2 where 

<f> = arctan ^- . In the real space notation the interac- 
tion term may be written as /fj = y^ i A^A^ where 
A, = X),5 ff(^) c i,i c i+<5-T- 9($) i s t ne Fourier transform of 
H(k) given as g(S) = jj-^2f:e lk ' s H(k). For example, for 
the following pair structure, 

r 1/2 if s = (±i,o), 

g(S) = i -1/2 if 5 = (0,±1), 
[ if otherwise, 

it leads to II(fc) = cosfc^ — cosk y . Obviously II(fc) = 
1.5264 cos20 requires in g(S) further neighbors beyond 
the four nearest neighbors. Here a numerical factor of 
1.5264 is introduced to normalize g(5) with J^s 5 2 (^) = 1- 
In terms of A^ the d-wave pairing susceptibility is defined 
as 

Xpp (i,r) = (T T A i (r)A+) , (2) 
where T T is the imaginary time ordering operator. 



It is well-known that for many strongly correlated elec- 
tronic models there is no obvious small parameter with 
which a systematic perturbative approximation can be 
made. In such a situation, a non-perturbative approach 
may be a good alternative for at least qualitative under- 
standing of the physics in those strongly correlated model 
Hamiltonians. In this paper we apply a non-perturbative 
sum rule approach to the model Hamiltonian (|l|), which 
is still not a fully controlled approximation in nature. 
Because of an extended nature of A^ for a d-wave pair 
and the increasing difficulty of the derivation in the real 
space representation, here we use the analogy from our 
previous studies for the repulsive |4(]]37[] and attractive 
Hubbard models. |36|] These studies show the important 
many-body modification with respect to the standard 
RPA and T-matrix approaches. The modification comes 
in two different places. First, the paring susceptibility 
is calculated through vertex function U pp instead of bare 
interaction strength U , which is constant in our approx- 
imation. Second, the pairing fluctuations are coupled 
to electrons, leading to the self-energy with UU PP form 
instead of U 2 . In particular, the latter structure is in 
agreement with the fact that there is no Migdal theorem 
for this problem, contrary to the case of electron-phonon 
interactions. And the renormalized interaction constant 
U pp is determined by the exact sum rule for the pairing 
susceptibility. By using the sum rule and the renormal- 
ized constant U pp , we determine in effect the Ginzburg- 
Landau parameters due to mode-mode coupling. Thus 
this approach is similar to the one-loop renormalization 
group approach within the Gaussian approximation and 
as a result the Mermin- Wager theorem is formally sat- 
isfied in two dimensions. EJ] This sum rule approach 
was systematically compared with the QMC simulations 
for the repulsive and attractive Hubbard models and the 
agreement was in a quantitative level both in one- and 
two-particle functions. We believe that this generic fea- 
ture of the many-body modification should also carry 
over to a Hamiltonian with d-wave pairing symmetry. 
In this modification, the pairing susceptibility and the 
self-energy can be written in Fourier space in terms of 
renormalized interaction strength V pp 



X PP (l) 



i + v PP x° PP (q) ' 

T 



x(k) = ~vv pp n 2 (k)-J2xpp(q)G°(q-k) , (3) 
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where the irreducible susceptibility is defined as 

x PP (q) 



^H 2 (fc)G°( (Z -fc)G (fc) 

k 



(4) 



The above expressions are reduced to the standard T- 
matrix approximation when V pp is replaced by bare V. 
We determine this renormalized constant by employing 
the exact sum rule for the d-wave pairing susceptibility 
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q 

= £s(%(<n<4 a Q, iC + CT ). (5) 

Note that the convergence factor is necessary in the sum 
rule, because away from half-filling Xpp{^ v m) decays like 
\jv m at large frequencies. In the previous studies for 
the repulsive and attractive Hubbard models, the right- 
hand side of the sum rule evaluated in the SDW and 
BCS mean-field ground states, respectively, was found 
to be in excellent agreement with the QMC values in 
the intermediate to strong coupling regimes. 0] In this 
paper the right-hand side of Eq. (Q) is also evaluated in 
the d-wave BCS mean-field ground state 

k k 

+ [^Y.m)u{k)v{k)] 2 , (6) 

fe 

where 



"fe - 2 ^ 1 E % ' ' 
fe 2 V 

% = \/(^-M) 2 + A 2 (fc). (7) 

Here A(fc) is the d-wave BCS mean-field gap. The chem- 
ical potential /x and the gap A(fc) are determined self- 
consistently through the number and gap equations for 
given V, T and n. 

Before starting we comment on some differences as- 
sociated with n(fc) = cos k x — cos k y and with H(k) = 
1.5264 cos 2<j> structure. The interaction Hamiltonian 
with the first structure depends not only on the angle 
of k but also on its magnitude. Thus the pairing interac- 
tion is always strongest at k, k' = (±7r, 0) or (0, ±7r), al- 
though the noninteracting Fermi surface can be far from 
these points. As a result two features occur for a particle 
density far away from half-filling (Detailed calculations 
were performed but not shown in this paper). First, 
the peaks associated with the precursor of the Bogoli- 
ubov quasiparticlcs occur asymmetrically with respect 
to the Fermi energy. Second, the locus of k points sat- 
isfying u> — e? + [i — ReY,(k,oj) = at to = can be 
substantially different from the noninteracting Fermi sur- 
face, thus strongly violating the Luttinger's theorem. For 
n(fc) = 1.5264 cos 2(f) which depends only on the angle, 
however, the above features disappear and the locus of 
k points satisfying the above equation is almost identical 



with the noninteracting Fermi surface. Near half-filling 
the differences are negligible. Throughout the calcula- 
tions, lattice spacing, h, and ks are set to be unity, and 
all energies are measured in unit of t. We used a dis- 
crete lattice as large as 128 x 128 in momentum space 
and performed the calculations by means of fast Fourier 
transforms (FFT) . Equations (||) and (^) are computed in 
terms of Matsubara frequencies and the analytic contin- 
uation from Matsubara to real frequencies are made via 
Pade approximants. p3fl In order to detect any spurious 
features associated with numerical analytical continua- 
tion, we performed real frequency calculations in paral- 
lel. Except for some spiky features in the real frequency 
formulation coming from the Lorentzian approximation 
of the non-interacting Green's function, the two results 
are almost identical. 

We begin by presenting the spectral functions along 
(0,0) - (vr,0) direction for V = -4t, n = 0.5 and 
T = 0.15t in Fig. [I]. Throughout the paper the den- 
sity is fixed at n = 0.5. The results for other densities 
are similar. For this parameter, V pp satisfying the sum 
rule is found to be — 1.87t significantly different from the 
bare value (V = —At). This shows the importance of the 
mode-mode coupling effect already in the intermediate 
coupling regime. Let us label wave vectors by (to7t/8, 0). 
Below the Fermi wave vector (to = 5), the main peak 
stays below the Fermi energy and the secondary peak 
grows in strength as k approaches kp . At the Fermi wave 
vector two peaks reminiscent of the Bogoliubov quasipar- 
ticles appear almost symmetrically with respect to the 
Fermi energy. Since electrons are still in the normal state, 
this is the precursor of the superconducting gap, namely, 
a normal state pseudogap. Above the Fermi wave vector, 
the main peak stays above the Fermi energy and the sec- 
ondary peak becomes stronger as k approaches kF- This 
result should be contrasted with that by Engelbrecht 
et al. [[mJ In the self-consistent T-matrix approximation 
with n(fc) = cosk x — cos k y form, these authors argued 
that along (0, 0) — (it, 0) direction the dominant peak of 
A(k, lo) never crosses the Fermi energy and bounces back 
towards the negative frequency. Their overall finding, 
however, is qualitatively different from our results. In the 
present calculations, the dominant peak of A(k, ui) even- 
tually passes through the Fermi energy for k far above the 
Fermi wave vector. At the noninteracting Fermi surface 
our spectral weight is strongly suppressed at the Fermi 
energy to become a local minimum, while the local min- 
imum of A(k,Lu) in their study is located significantly 
away from the Fermi energy. Presumably this is due to 
the features associated with n(fc) = cos k x — cos k y struc- 
ture mentioned in the previous paragraph and also due 
to the self-consistent approximation that does not take 
the vertex function and the Green's function at the same 
level of approximation. {3 

In Fig. ^ the spectral function is shown along the Fermi 
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wave vectors at low temperature (T = 0.15i). The angle 
<j> is denned as a,rcta,n(k y / k x ) along the noninteracting 
Fermi surface. At <j> = 0° the magnitude of the pseu- 
dogap is largest and as 4> increases it progressively de- 
creases, leading to strong momentum or angle depen- 
dence in the size of the pseudogap. Since along the diag- 
onal directions the pairing interaction vanishes, the pseu- 
dogap completely closes and a sharp quasiparticle peak 
appears at (j> = 45°. In Fig. || the magnitude of the 
pseudogap (circles) is plotted as a function of angle (j> 
for T — 0.1 ht along with the ground state gap symmetry 
A(</> = 0°)cos2</> (dashed curve). The angle dependence 
of the normal state pseudogap is well fitted by d-wave 
symmetry, consistent with ARPES experiments for un- 
derdoped cuprates. 

At higher temperatures a drastic change is found in the 
spectral function. For T = 0.225i the pseudogap closes 
well below 45°, as shown in Fig. [|. The local minimum 
at the Fermi energy disappears at <j> = 31.0° (26.6° is just 
on the crossover). Beyond this angle, the quasiparticle- 
like peak appears and thereby the Fermi surface is par- 
tially restored. We define this angle as critical angle <p c . 
In Fig. |^ <j> c is plotted for different temperatures (open 
circles). Below T = 0.175i the pseudogap is found every- 
where along the Fermi wave vectors except at <fi = 45°, 
thus the Fermi surface is destroyed everywhere except 
along the diagonal directions. With increasing temper- 
ature (up to T = 0.29t), however, <\> c becomes smaller 
than 45°. Thus, the pseudogap region shrinks and at the 
same time the Fermi surface grows from <j> = 45° up to 
the critical angle. Above T = 0.29i the whole Fermi sur- 
face is completely restored in spite of some broadening 
of the spectral function due to interaction as well as fi- 
nite temperature. This feature is qualitatively consistent 
with ARPES for underdoped cuprates. [|| We can theo- 
retically calculate the temperature dependence of <p c . By 
using the Ornstein-Zernike form of the pairing correlation 
function and taking the classical fluctuations (iq n = 0), 
pEf the scattering rate at the Fermi energy is found to be 
proportional to II 2 (fc)^/(T, where £ and £t are pairing 
correlation length and thermal de Broglie wave length de- 
fined as £t = Vf{k)/T, respectively. For the pseudogap 
to disappear, the scattering rates should be much smaller 
than unity II 2 (£;)£/£t <C 1, allowing us to define a critical 
angle £ = vf{4>c) /T /H 2 ((j) c ). In Fig. || this critical angle 
is also shown as stars. For the best fit near 0° angle, a 
numerical factor of 1.06 multiplies £. Although a small 
deviation is found near 45°, the overall magnitude and 
shape are in reasonable agreement throughout the whole 
angle. For a d-wave model, the condition for the appear- 
ance of a pseudogap in a given momentum depends not 
only on the anisotropy of the Fermi velocity (which is the 
only relevant condition in the attractive Hubbard model) 
but also more importantly on an angle dependent form 
factor n 2 (fc). 



In Fig. |6j we show the spectral function at the Fermi 
wave vector and the density of states for different tem- 
peratures. As the temperature is increased, the spectral 
weight starts to fill in the pseudogap and at T = 0.3i the 
precursor of the superconducting gap completely disap- 
pears as shown in Fig. ^(a). At this temperature, how- 
ever, the pseudogap still persists in the density of states 
and at higher temperature (T = 0.4t) it finally disappears 
as shown in Fig. |](b). Except close to half- filling, the 
pseudogap appears at higher temperature in the density 
of states than in the spectral function. Compared with 
the s-wave case, the density of states is suppressed lin- 
early near the Fermi energy, a reminiscence of N(uo) ~ u) 
in the superconducting state. Both in the spectral func- 
tion and the density of states, the pseudogap disappears 
in a manner that the spectral weight fills in the pseu- 
dogap instead of closing it with increasing temperature. 
This feature is also consistent with ARPES for under- 
doped cuprates. This may suggest that in our approach 
phase fluctuations (spin-wave type) rather than ampli- 
tude fluctuations are mainly responsible for the pseudo- 
gap formation, although the present approach includes 
both. Like in the s-wave case, a pseudogap also appears 
in the density of states when the characteristic pairing 
frequency v c is equal to or smaller than temperature. 
This corroborates the origin of the d-wave pseudogap, 
namely, growing d-wave paring fluctuations in the low 
temperature renormalized classical regime of the low di- 
mensional system. The calculated T* follows the same 
trend as the mean-field critical temperature Tmf- (To 
be more precise, T* is approximately half of Tmf for 
most of the densities.) As noted in Ref. [26| ], near a 
point with high order parameter symmetry (half-filling 
in that paper), the transition temperature T c decreases 
while the pseudogap temperature increases along with 
Tmf- As a result it leads to a large pseudogap regime, 
consistent with the phase diagram in the underdoped side 
of cuprates. 

There are several advantages in our formulation. First, 
the pairing fluctuation sum rule Eq. (^|) is exactly sat- 
isfied (by construction). Through this sum rule, the 
Mermin- Wagner theorem is formally fulfilled and the 
strength of pairing fluctuations is properly constrained 
within the Gaussian approximation. This latter fea- 
ture is crucial in our formulation, because an approxi- 
mate treatment of pairing fluctuations without constrain- 
ing the strength can easily overestimate or underesti- 
mate the magnitude of fluctuations particularly in low 
dimensional systems. The Mermin- Wagner theorem is 
also satisfied in some other approaches such as the self- 
consistent T-matrix (FLEX) and the "pairing" approx- 
imation schemes. Second, there is an exact relation 
between one-particle (self-energy, Green's function) and 
two-particle (interaction term) functions: 
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lim Y, ^$,iu n )G{k,iu n )e-^ T = V(A+A,} , (8) 

r^0~ IS _f — ' 

where here the self-energy includes the Hartree-Fock 
term. When the interacting Green's function G(k, iuj n ) is 
replaced by the noninteracting one G°(k, iuj n ) in Eq. (|^), 
it is exactly satisfied. With G(k,ioj n ), the difference be- 
tween the left- and right-hand sides is less than 6% for 
all temperatures studied. In this paper the Kosterlitz- 
Thouless phase transition |l6| and its fluctuation effect 
near Tkt = T c have been neglected, since the present for- 
mulation is inadequate to describe the topological nature 
of the phase in two dimensions and vortex anti-vortex 
binding- unbinding physics. Since for similar parameters 
for the attractive Hubbard model the QMC results indi- 
cate that Tkt ~ 0.05i j23] three times smaller than the 
lowest temperature in our calculations, we do not expect 
in the present results any significant influence from vor- 
tex phase fluctuations. Our approach is valid for weak to 
intermediate coupling and for temperature not too deep 
in the pseudogap regime. At very low temperature, even 
the Ginzburg-Landau functional form itself may change, 
for instance, a possible crossover of the dynamical criti- 
cal exponent from z=2 to other value, jlTj and eventually 
vortex phase fluctuations may come into play. 

In summary, we have studied spectral properties and 
pseudogaps in a model with d-wave pairing symmetry 
by using a non-perturbative sum rule approach. We ap- 
plied to this model our previous experience of many- 
body theory in the repulsive and attractive Hubbard 
models. The magnitude of the normal state pseudogap 
shows strong angle dependence well fitted by cos 2<p form 
at low temperature. With increasing temperature the 
pseudogap closes at some critical angle <f> c and beyond 
this angle a single quasiparticle-like peak appears. The 
resulting Fermi surface is strongly temperature depen- 
dent. With increasing temperature, the pseudogap re- 
gion shrinks and at the same time the Fermi arc grows 
from <j> = 45° to <fi c . Both in the spectral function and 
density of states the pseudogap disappears in a manner 
that the spectral weight fills in the pseudogap instead of 
closing it with increasing temperature. All these features 
are qualitatively consistent with ARPES for underdoped 
cuprates. The pseudogap is caused by growing d-wave 
critical pairing fluctuations in the low-temperature clas- 
sical renormalized regime of the low-dimensional system, 
as in the repulsive and attractive Hubbard models. We 
argue that although the real critical behaviors and crit- 
ical exponents are governed by the vortex phase fluc- 
tuations close to the Tkt, the initial growth of pairing 
fluctuations can be driven by spin-wave phase fluctua- 
tions, leading to the normal state pseudogap formation. 
The behavior of the spin-wave type phase fluctuations 
belonging to the 0(n = oo) universality class |58| can 
be qualitatively different in two dimensions (particularly 
near a point with high order parameter symmetry) from 



that of the 0(2) Kosterlitz-Thouless vortex phase fluctu- 
ations. 
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FIG. 1. Spectral function along the line (0,0) — (n, 0) for 
V = —4, n = 0.5 and T = 0.15. The figures for m = 3 - 8 
are shifted vertically by 0.25. m = 5 corresponds to the non- 
interacting Fermi surface. 



FIG. 3. The magnitude of the pseudogap for V = —4, 
n = 0.5 and T — 0.150 (open circles). The angle <f> is defined 
as arctan(fe I/ /fca;) along the noninteracting Fermi surface. The 
dashed curve is A(<f> = 0°) cos 20. 
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FIG. 2. Spectral function along the noninteracting Fermi 
surface for V = —4, n = 0.5 and T = 0.15. The angle <j> 
is defined as arctan(fcj / /fc a: ) along the noninteracting Fermi 
surface. The figures for cf> > 0° are shifted vertically and 
horizontally both by 0.5. 
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FIG. 4. Spectral function along the noninteracting Fermi 
surface for V = -4, n = 0.5 and T = 0.225. The angle (f> 
is defined as a,Tctan(k y / k x ) along the noninteracting Fermi 
surface. The figures for <f> > 15.5° are shifted vertically and 
horizontally by 0.25 and 0.5, respectively. 
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FIG. 5. Critical angle <f> c as a function of temperature for 
V = —4, n = 0.5 (open circles). The stars are the critical 
angle estimated by the expression given in the text. 




FIG. 6. (a) The spectral function at the noninteracting 
Fermi surface and (b) the density of states for V = —4 and 
n = 0.5 at different temperatures. The solid, dotted, dashed, 
long-dashed, and dot-dashed curves correspond to tempera- 
tures from top to bottom. 
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